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5? : ABSTRACT 

Combining heavy quark effective tfieory and tfie cliiral Lagrangian ap- 
proacli we investigate radiative decays of pseudoscalar D mesons. We first 
reanalyse D* -^ D'y decays within the effective Lagrangian approach using 
heavy quark spin symmetry, chiral symmetry Lagrangian, but including also 
the light vector mesons. We then investigate D — > V'j decays and calculate 
the D^ -^ K*^'-f and D^^ -^ p'^'-f partial widths and branching ratios. 



1 .Introduction 

The list of D meson decay rates is rather long and further study of their 
decays would eventually help to better understand their features. There has 
not yet been any experimental evidence for radiative decays of D mesons, 
while D* radiative decays are known to be important. The D* decays [|I], 
can be described in a model independent framework which incorporates 
the appropriate constraints on the decay amplitudes. The combination of 
heavy quark effective theory and chiral Lagrangians have been extensively 
studied and applied to many D mesons decays [^, ^, |^, |^, 0, |], ^, |TD|, |TT], 
|1^, 0. Wise [1^ has proposed an effective Lagrangian to describe, at low 



momentum, the interactions of a meson containing a heavy quark with the 
light pseudoscalar mesons tt, K, rj. Two kinds of symmetries characterize 
the effective Lagrangian: the heavy quark SU{2) spin symmetry and the 
non-linearly realized SU{3) ® SU{3) chiral symmetry in the light sector, 
corresponding to spontaneous symmetry breaking of the chiral group to the 
diagonal SU{3)v- Due to the rather large masses of the D mesons, the 
inclusion of resonances with masses below the D mesons seems necesssary 



n|,[T|,[Tl 



In this paper, following the requirements of heavy quark and chiral sym- 
metry, we develop a framework for the description of heavy and light pseu- 
doscalar and vector mesons. In section 2 we write down the most general 
Lagrangian in the limit of exact heavy quark and chiral symmetries. Section 
3 is devoted to the higher order odd parity Lagrangian, which also describes 
the decay D* -^ D^ , and we reinvestigate this decay in order to learn more 
about the couplings in the chiral Lagrangian. In section 4 we analyse the 
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weak Lagrangian. Finally, as an example of the use of our model, we calculate 
the D —^ V'-f radiative decays in section 5. 

2. The chiral Lagrangian technique and heavy quark Umit 

The strong interaction meson Lagrangian for the light pseudoscalar octet 
and heavy pseudoscalar and vector triplets in the chiral and heavy quark 
limits was first written down by Wise |T^ (see also 0]). The electromagnetic 



interactions between these mesons was described in [|T], ^|^. The octet of light 
vector mesons was included in the Wise Lagrangian |I^ later by Casalbuoni, 
et al. |T^ as the gauge particles associated with the hidden symmetry group 
SU(3)/^ |T^. The next step is to provide a common description of both 
the light and heavy pseudoscalar mesons, which also includes both light and 
heavy vector mesons and the electromagnetic interactions. In this section 
we present the strong and electromagnetic Lagrangian for the description of 
both light and heavy pseudoscalar and vector mesons. 
The light pseudoscalar mesons are described by the 3x3 unitary matrix 



in 
u = exp(— ) (1) 

where / ~ 132 MeV is the pion pseudoscalar pion constant and 11 is the 
pseudoscalar meson unitary matrix defined as 



n 



v/2 ^ v/6 


TT+ 


K+ \ 


7T~ 


V2 ' Ve 


K' 


\ K- 


K' 


>/ 



(2) 



The octet of light vector mesons is described by the 3x3 unitary matrix 



P/^ 



■ 9v 



(3) 



where gv (— 5.8w2/a with a = 2 in the case of exact vector dominance) 
is the couphng constant of the vector meson self-interaction |T^ and p^ the 
vector meson unitary matrix 



/ 



Pf^ 



V2 



P 






v/2 






(4) 



In the following we will also use the gauge field tensor F^y{p), defined as 



Ft,u{p) = d^p^ - d^p^ + [p^, p^] 



(5) 



The heavy mesons are Qq'^ ground states, where Q is a c (or b) quark and 
q^ = u, q"^ = d and q^ = s. In the heavy quark limit they are described by 
4x4 matrix H^ (a = 1, 2, 3) 



H. 



l + ^){P*r - Pal,) 



(6) 



where P* and Pa annihilate, respectively, a spin-one and spin-zero meson 
Q(f of velocity v^. The creation operators P^^I ^^^ P^ occur in |jl4 



Ha 



H-t.,0 



I'Hli 



(^:i7'^+^i75)^(i+? 



1 

■2' 



(7) 



Following the analogy in refs. [^, |T^ we introduce two currents: 



V^ = ]^{u^D^u + uD^u^) (8) 

and 



A^ = -{u^D^u-uD^u^) (9) 

where the covariant derivatives of u and u^ are defined as 



D^u = {d^ + B^)u (10) 



and 



D^u^ = {d^ + B^)u^ (11) 



while 



B^ = leB^Q (12) 







g = I -I I (13) 



3 , 

-I 



and B^ is the photon field. To insure that the vertices D^'^D^^ and D*^''D*^'y, 
or those with D replaced by B in the case of the b quark, are absent, we define 
the covariant derivative for the heavy meson field as 



D^Ha = id^ + V^-teQ'B^)Ha (14) 

with Q' = 2/3 for c quark (—1/3 for b quark). With these definitions we can 
finally write down the even parity strong and electromagnetic Lagrangian for 
heavy and light pseudoscalar and vector mesons: 



+ ^gTr[H,-f,-f,{AnbaHa]+^f3Tr[H^,v^{V^ - p^),aHa] (15) 

2f^a 
with 



f2 

Ciight = -y{tr(^^^^) + atr[(V^-p^)2]} 

+ Atr[F,,ip)F''''ip)] (16) 

This Lagrangian is invariant under the following gauge transformation: 

u -^ go{x)ugl{x) 
u^ -^ go{x)u'< gl{x) 
V^ -^ 9o{x)Vf,gl{x) + go{x)df,gl{x) (17) 



Pti -^ 9o{x)pf,go{x) + go{x)d^gl{x) 
Bf, -^ go{x)Bf,gl{x) + go{x)df,gl{x) 

where go{x) = exp{ieQX{x)). The last transformation (|l^) together with 
([T2D -(|l3|) imply, of course, the usual gauge transformation for the photon 
field: 

B^ ^ B^-d^Xix) (18) 

In equation (p!5D g and (3 are constants which should be determined from 
experimental data |jl|, 0, |Tl|, |1^, |13|. The constant a in (p!5D-(p!6D is in principle 
a free parameter, but we shall fix it by assuming exact vector dominance [|l^ 



for which a = 2. With exact vector dominance there are no direct vertices 
between the photon and two pseudoscalar mesons, so that the pseudoscalars 
interact with the photon only through vector mesons. 

The electromagnetic field can couple to the mesons also through the anoma- 
lous interaction; i.e., through the odd parity Lagrangian. Even with the 
PP'y direct vertices absent in Cught due to the choice a = 2, direct PV'y 
vertices are present in the odd parity Lagrangian. We write down the two 
contributions which are significant for our calculation: 

C'-l = -^^e^-'Tr{d,p^d^ppU) (19) 



4?. = -4ev^^e^-^Tr({a,p.,n}Q9.i?^) (20) 



6 



Equation ([T9|), together with vector dominance couphngs 

Ly., = -m'y^B,{p'^ + \u^-^<^n (21) 

gv 6 6 

which come from the second term in (p!6D, describe the electromagnetic in- 
teraction assuming vector-meson dominance, while the direct photon-light 
vector meson-pseudoscalar interactions are contained in (^). The contri- 
butions to the odd Lagrangian (^) and ( PPD arise from Lagrangians of the 
Wess-Zumino-Witten kind |16|, ^. 

In the rriq —>■ and ttlq -^ oo limit {rriq and rriQ are the masses of the light 
and heavy quarks, respectively) the strong and electomagnetic interactions 
of heavy and light pseudoscalar and vector mesons are thus described by the 
even Lagrangian (|l5D-(0) and by the odd Lagrangian ([T9|)-(pOD. However 



the D*D'~f vertices are not included in the above Lagrangian since it is of the 
anomalous type. The terms responsible for it are of higher order I/tuq and 
they will be introduced in the next section. 

3. Higher order odd Lagrangian for heavy mesons 

In our approach vector meson dominance describes the couplings of light 
quarks and photons through the higher dimensional invariant operator 

C, = iXTr[HaCX,,F^'''{p)abHb] (22) 

In this term the interactions of light vector mesons, heavy pseudoscalars or 
heavy vector D mesons are also present. The light vector meson can then 
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couple to the photon by the standard vector-meson dominance prescription 
(|2T|). These terms effectively describe the light quark- photon interaction 
inside the charmed (or beauty) mesons. 

The coupling A can be independently determined either from D* decays into 
D-n and D7 [|T], ^ ^ |] , some ratios of which have been measured [|T8|, |19| or 
from semileptonic decays of D mesons |Tl], |12|, |13] . 



The heavy quark-photon interaction is genereted by the term 



£2 = ->:eTr[Ha(y^,F^^{B)Ha] (23) 

According to quark models the parameter | A' | can be approximately related 
to the charm quark magnetic moment via l/(6r7ic) |l|, ^, ^ |^. In order to 
reduce the error in determinig the couplings we shall reanalyze the decays 
D* -^ D-K and D* -^ D'-f. Experimentally one measures [|I^] the branching 
fractions i?° = r(L)*° -^ D^-f)/T{D*^ -^ DV) = 0.572 ±0.057 ±0.081. and 
/?+ = T{D*+ -^ D+-i)/V{D*+ -^ D+7r°) = 0.035 ± 0.047 ± 0.052. Using our 
Lagrangian these branching ratios are 



9 Sg Pn 



i?; = M7rfa^M{- + --yC-ir (24) 



^ 9 Sg Pt 



i?+ = 6A7^faEMi-~=:-n^r (25) 



To determine \/g and X' /g, the square-roots of the left-hand-sides of eq. 
(|24D-(p5|) have to be taken, which introduces an ambiguity in the resulting 



coupling constants, evaluated. The experimental errors in the branching 
fractions are somewhat large but the masses of the particles involved are 
relatively well known. We have used the standard formulae flSl 



/ = f{x) + ^cr'f"{x) (26) 



aj = a\f{x)\ (27) 

for /(x) = y/x (x = R^ or R'i^), where x, f and a, af are the mean values 
and standard deviations of the respective distributions. Using expressions 
(PBD-(P^) instead of their linearized versions is important due to the large 
experimental error in i?;|. From that data we then find 



I- --I 
9 ^9 



(0.863 ±0.075)Ger"^ (28) 



and 



|--^-| = (0.089 ± 0.178)^6^-^ (29) 

9 '^9 

The two errors in eq.(p8[)- ([29| ) can in principle be correlated. The main 
sources of correlations are the experimental efficiencies in the detection of 
the 7r° and 7. Since the contributions of the efficiencies to the net errors 
are small ||IP[, possible correlations were neglected and the errors in the de- 



termination of single \/ g and \' / g were combined in quadrature. There are 
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four solutions for the ratios Xjg and \' / g: a) \/\g\ = /i;(0.77 ± 0.19)GeV ^, 
X'/\g\ = k{0.35 ± 0.12)GeV-i and b) X/\g\ = k{0.95 ± 0.19)GeV-\ X'/\g\ = 
A;(0.23±0.12)GeV~^, where k can be +1 or —1. For the individual determina- 
tion of A, A' and g, one has to fix one parameter. Choosing |A| = (0.60±0.11) 
GeV~^ as determined in |T^ minimized the errors in these quantities. We 
present in Table |I] the four possible solutions of (p8D-(p^) for A and \g\, to- 
gether with the combinations A' + 2A/3 and A' — A/3, which will be used in 
our calculations in section 5. 





A [GeV-i] 


A' [GcV-i] 


l<?l 


(A' + f A) [GeV-i] 


(A' - lA) [GeV-i] 


1 

2 
3 

4 


+0.60 + 0.11 
-0.60 + 0.11 
+0.60 + 0.11 
-0.60+0.11 


+0.29+0.12 
-0.29 + 0.12 
+0.15 + 0.09 
-0.15 + 0.09 


0.82 + 0.20 
0.82 + 0.20 
0.66 + 0.13 
0.66 + 0.13 


+0.71 + 0.18 
-0.71 + 0.18 
+0.57+0.12 
-0.57 + 0.12 


+0.07 + 0.15 
-0.07 + 0.15 
-0.06 + 0.12 
+0.06 + 0.12 



Table 1: Four possible solutions from (28)-(29) with A determined by [13]. 



The experimental value l^^l = 0.57±0.13 (see for example \I^\ and references 
therein for a discussion of these experimental values) and the approximate 
validity of the equation |A'| ~ l/(6mc) with rric — 1.5 GeV slightly favour 
solutions 3) and 4). 

Our approach is different from |]T]|, ||^ and |^, since we do not use any quark 
model prediction for the parameter A' but treat it on an equal footing with 
the parameter A, so that both are considered as purely phenomenological. 
Nevertheless we were able to obtain reasonably good precision in the deter- 
mination of model parameters. 
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4. Weak Lagrangian for light and heavy mesons 

In addition to strong and electromagnetic interactions, we must also address 
the weak decays within these scheme. We will follow the approach of [O, P 



and use an effective current between the heavy mesons and the light mesons. 
The weak current is Lq^ = ga7^(l —l5)Q and it transforms as (3^, 1r). The 
lowest dimension operator with the same transformation properties but with 
meson degrees of freedom is 

jQ'a = l^aTr[r{l-l,)H,ul] (30) 

+ aiTrl^.Hhip" - V^tcul] + ■■■ 

where the ellipsis denote terms vanishing in the limit niq —* 0, rriQ — » oo or 
terms with derivatives. 
The light mesons decay constants fpy are defined by the usual relations 



\ab 



< o\j,l\omip) > = ^fp^P, 



\f 



<0|<(0)|V^,(e,p)> = fv^mye, (31) 

Here A"*, i = 1, . . . , 8, a,b = 1,2,3 are the usual eight Gell-Mann 3x3 
matrices, normalized as Tr{XiXj) = 26ij. In the chiral limit rriq —>■ the 
above decay constants are related through fp = jy j \fa = f for all P, V. 



Similarly we can define the heavy meson decay constants by [14 



<0|Jq;(0)|L>^(p)> = -tfoS^'mnv^ 
<0\jQliO)\D*\e,p)> = t/n^d'^'niD^e, (32) 
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In the heavy quark hmit tuq ^ oo we have mD = rriD* -^ oo and Jd = fo* ■ 
The constant a in eq. (^) can be fixed in this hmit by taking the matrix 
elements of Jj^ between the heavy meson state and the vacuum, with the 
result 



a = fp^/rrip (33) 

Unfortunately, up to now, there does not exist either a theoretical prediction 
or experimental data for the other parameter, ai, in the current (|30|) . 

The situation is different and better in the light sector, where a well known 
prescription exists, for deriving the weak current directly from the strong 
Lagrangian [^. Since the quark weak current Lg^^^ = gfe7^(l — 75)(3'a must 
transform as {3l, 3r) the light meson weak current with these transformation 
properties can be obtained from the Lagrangian (|TBp. Of course, one has to 
properly define the covariant derivative for pseudoscalars, enlarging the gauge 
group to include the VT-boson contributions [^]. The resulting light meson 
part of the weak current is 

J/ = ifu[A'' + a{V''-p^)W (34) 

The part of the weak Lagrangian for the pseudoscalar and vector, light and 



heavy mesons, which we will use, can be written as p2|, ^, |^ 



C%^{Ac = As = l) = -^VudV:A ai{ud)'i,_^{sc)v-A,^.+ 

a2{sd)'^^^{uc)v-A,fj. ] (35) 
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where Vud, etc. are the relevant CKM mixing paremeters, while ai and 02 
are the QCD Wilson coefficients, which depend on a scale yU. One expects the 
scale to be the heavy quark mass and we take yU ~ 1.5 GeV which gives ai = 
1.2 and 02 = —.5, with an approximate 20% error. In the factorization model 
the quark currents are approximated by the corresponding meson currents 
defined in eqs. (|30| ) and (pi]): 



{qaQry_^ ^ g,7'^(l-7')Q^jQ^ (36) 



aj'v-A = qi>Yii-i')qa^V,, (37) 



Many heavy meson weak nonleptonic amplitudes p2|, |23|, ^ ^ have been 
calculated using the factorization approximation. It has been shown in [p^] , 
however, that for some of the D meson decays there are rather important 
final state interactions and the factorization approximation can be improved 
by the inclusion of the SU{3) symmetry breaking effects ||26[ . 



The authors of ref. [^ have classified the weak nonleptonic decays into 
three classes: decays determined by ai only (class I), decays determined by 
02 only (class II) and decays where ai and 02 amplitudes interfere (class III). 
Factorization can therefore be tested in several ways. There are the following 
two cathegories of decays: " quark decays" , in which the heavy quark decays 
while remaining antiquark acts as a spectator, and "annihilation processes" 
in which heavy and light quarks annihilate and two new quarks are created. 
For the annihilation processes the factorization approximation is usually only 



a small contribution [23, 24, 
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We are forced to use this approximation in our calculations, since there are 
no better approaches developed so far for nonleptonic weak decays. 

5. -D — > V'j decays 

The simplest radiative decays of D mesons are into a light meson and a 
photon. Since the process D — > P7 (P is a light pseudoscalar) is forbidden 
due to the requirement of gauge invariance and chiral symmetry [|2^ , as well 
as angular momentum conservation, we will concentrate on the D —>■ Vj 
(V is a light vector meson) decays. Although there are not yet any data on 
such processes we can predict the partial widths and branching ratios. We 
consider the only two processes which are possible at tree level and are not 
Cabibbo supressed, namely D^ -^ K*^'j and D^~^ -^ P^l- Both processes 
have contributions from the odd-parity interaction Lagrangian. The second 
one has, in addition, a direct emission term, due to the charged initial and 
final mesons. 

Regarding the anomalous term, there are two contributions which are im- 
portant. The photon can first be emitted from the D^ (-D*^) meson, which 
becomes a D*^ {D*"^) and then D*^ (/}*''+) decays weakly into A"*o (p+). 
The other contribution comes from the weak decay of D^ {D^^) first into an 
off-shell K^ (tt^), which then decays into K*^'^ {p^l)- Both contributions are 
proportional to 02 (o-i) (pS]). For the description of this amplitude we need 
the D*D'y and K*K'y (pTC'y) couplings and these couplings were obtained in 
the previous section. 

The amplitude for the D^ -^ K*^'y is 
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A{D%pn.) ^ K*\pK^.)^{q)) = e^KdV;>2 (38) 



where 



.(1) _ rr^ ^ , r^ \ fDofKomjjo Sy^moo 



^ 9v {mjjo - mj^o) 3/ 

+ 4(A + -A)/d*o//^*o— ^ 2 — ^V^DorriD'O (39) 

Similarly, the amplitude for D^^ -^ p+7 is 



+?r(^^ m (f* f* - ^xMkl]] 

Pp-H 



where 



^(1) ir^ ^ ^n \ fD-U^h 4V2mD» 






+ 4(Y--A)/,,»/ - ' V^^^E^^^^ (41) 



and 



Cdv7 - /c^/p (42) 
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In our numerical calculations we used the following numerical values Cvvn = 
.423, Cvwy = -3,26.10-2, p8|, |2|, gv = 5.8 ^, f - f^ = 132 MeV, and 
the other decay constants fpy were taken from [^. It is straightforward to 
calculate the decay widths. The result, of course, depends on which numerical 
value we take for (A' ± 2A/3) and (A' — A/3). The numerical predictions for 
the decay widths and branching ratios are given in Table |^, where all the 
four possible choices of paramaters (see Table |l|) are considered. 



Solution 


D' -^ K*^^ 


D'+ ^ p+7 


r lO-'^eV 


BR 10-4 


r lO-'^eV 


BR 10-4 


1 
2 
3 

4 


7.8 ±2.5 
0.8 ±0.7 

5.9 ±1.4 
0.3 ±0.3 


4.9 ±1.6 
0.5 ±0.4 
3.7 ±0.9 
0.2 ±0.2 


4.9 ±2.0 

7.6 ±3.7 
7.0 ±2.9 

4.7 ±1.7 


3.3 ±1.3 

5.1 ±2.5 
4.7 ±1.9 

3.2 ±1.1 



Table 2: Different possible predictions for the D^ — > K*^'j and -D'*"'" -^ 
p+7 decays. The errors are due to uncertanties in the determination of the 
combinations (A' ± 2A/3) and (A' — A/3). 



An interesting feature can be seen from Table ^: A not very precise mea- 
surement of the D^ —>■ K*^'~f decay rate is sufficient to differentiate between 
solutions l)-3) and solutions 2)-4), which are predicted to be of one order 
of magnitude different. Unfortunately, due to the much larger branching 
ratio for the weak decay D^ -^ K*^7t^ and the difficulty in differentiating 
the photon from the vr" in this energy range, the decay D^ —>■ K*'-f has not 
been seen by the ARGUS collaboration [Q. The situation is similar for the 
detection of D^^ -^ p^'-f by the ARGUS collaboration: due to the very small 
branching ratio, low detector's acceptance (37 events have to be measured) 
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and poor mass resolution, the ARGUS data are not likely to find this de- 
cay. But, hopefully, some experimental signals for these radiative decays will 
come from the CLEO data. The experimental measurement of this branch- 
ing ratio will determine the relative sign between the first and the second 
contributions, giving in such a way new information on the parameters A, A' 
and g. 

In conclusion, we used both chiral symmetry and heavy quark symmetry 
to obtain an effective strong, EM and weak Lagrangian for the description 
of both light and heavy pseudoscalar and vector mesons. In this framework 
we reanalyized the D* strong and radiative decays, obtaining without any 
reference to quark models, a good determination of some of the parameters 
in the effective Lagrangian. Within the same framework and with these 
values for the parameters we calculated the D — * V-^ decay widths, providing 
numerical predictions. These results can be used to test the validity of the 
approximations that were made in the context of the heavy quark effective 
theory. At the least, our numerical results are reasonable estimates and 
provide some guidance. In the framework developed here other D meson 
radiative non-leptonic decays {D — > FP7 or PP'y) can also be calculated [^], 
giving estimates for future experiments and further tests of the applicability 
of HQET. 

Acknowledgement. This work was supported in part by the Ministry of 
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German Ministerium fiir Forschung und Technologic under contract No. 06 
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Note added. After completing of this work the preprint of P. Jain, A. Momen 
and J. Schechter (|hep-ph/9406338| ) , where some of these results have been 



independently obtained, came to our attention. Their approach is quite 
similar to ours in section 3, but they fix the parameter A' to be proportional 
to l/rric, while we leave it free. Also we noticed a preprint by H.Y. Cheng, 
C.Y. Cheung, G.L. Lin, Y.C. Lin, T.M. Yan and H.L. Yu (|hep-ph/9407303| ). 



In this preprint the D —>■ K*'~f decay rate was estimated using an effective 
electromagnetic and weak Lagrangian developed from quark diagrams for 
bd -^ cu'j. They then made the replacement b ^ c and c — *> s. Their result 
is comparable with our second case in Table |^. 
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